Clustered Error Correction of Codeword-Stabilized Quantum Codes 



O 

(N 



oo 



43 

i 



> 

00 

m 
O 

o 

o> 
o 



- I - < 



Yunfan Li, 1 Ilya Dinner, 1 and Leonid P. Pryadko 2 

1 Department of Electrical Engineering, University of California, Riverside, California 92521, USA 
2 Department of Physics & Astronomy, University of California, Riverside, California 92521, USA 

(Dated: March 9, 2010) 

Codeword stabilized (CWS) codes are a general class of quantum codes that includes stabilizer 
codes and many families of non-additive codes with good parameters. For such a non-additive code 
correcting all t-qubit errors, we propose an algorithm that employs a single measurement to test all 
errors located on a given set of t qubits. Compared with exhaustive error screening, this reduces 
the total number of measurements required for error recovery by a factor of about 3 . 
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Quantum computation admits polynomial complexity 
for many classical algorithms believed to be hard |], ||. 
To preserve coherence, quantum computations must be 
protected by quantum error correcting codes ^] . Sta- 
bilizer codes 0] represent an important class of quan- 
tum codes that can be encoded and decoded in polyno- 
mial time. Recent Refs. ^, ^ introduce a larger class of 
codeword-stabilized (CWS) codes. It includes important 
code families, such as the stabilizer codes and generally 
non-additive union stabilizer (USt) codes ||. CWS codes 
have a broader range of code parameters which can be 
superior to those of any stabilizer code|^-^]. 

The most important advantage of the CWS codes is 
their close relation with the classical codes. In particular, 
a qubit CWS code Q can be mapped onto a classical bi- 
nary code C, with the quantum Pauli errors also mapped 
into some binary error patterns H. This way, within 
CWS framework, quantum code design can be reduced 
to classical codes and employ the wealth of different tech- 
niques developed for the latter. 

On the other hand, quantum error correction must pre- 
serve the original quantum state in all intermediate mea- 
surements, and therefore is more restrictive than many 
classical algorithms. Thus, design of CWS codes must 
be complemented by an efficient non-damaging quantum 
error correction algorithm. In this paper, our main goal 
is to address this important unresolved problem. 

We consider a general non-additive CWS code 
((n, K, d)) of distance d which encodes K quantum di- 
mensions into a if-dimensional subspace of the Hilbert 
space of n qubits. This code detects all errors that cor- 
rupt up to (d— 1) qubits, and corrects all errors corrupt- 
ing t = [(d — 1)/2J or fewer qubits. As a benchmark for 
our study, we consider generic algorithms that project 
a corrupted code state into different subspaces. This 
brute-force technique is similar to the exhaustive error 
screening in nonlinear classical codes, and requires up to 
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To reduce the number of such measurements, we first 
design an error detection algorithm for USt codes JsJ . In 
CWS framework, the classical code C associated with the 
USt code Q = ((n,m2 k ,d)) is decomposed as a group 
Co of 2 k codewords shifted by m binary "translation" 
vectors. We prove the following 

Theorem 1 For a USt code of length n, with a group 
of size 2 k and dimension K = m2 k , an error-detecting 



measurement requires no more than 2m(n 
two-qubit gates. 



k){n + 3) 



Then, for a general CWS code Q, we propose an error- 
correcting method that simultaneously screens all 4* dif- 
ferent errors located on any given subset of t qubits, by 
designing an auxiliary USt code which uses binary maps 
of these errors as generators of the group Co, and the 
codewords of the associated classical code C as transla- 
tions. This requires only (™J — 1 measurements to screen 
all groups. Once the corrupted qubits are located, we 
need up to 2t extra measurements to find the actual er- 
ror within the group. Overall, this reduces the number 
B(n,t) of measurements about 3* times to 



N(n,t) = ) +2t- 1. 



(2) 



Our main result is summarized as 



measurements to screen all errors of weight t or less. 



Theorem 2 Consider any t-error correcting CWS code 
of length n and dimension K . Then this code can correct 
errors using at most N(n, t) measurements, each of which 
requires at most 2K(n — l)(n + 3) two-qubit gates. 

Definitions. Throughout the paper, we use the 
Hilbert space %® n = H 2 to represent any n-qubit 
state. Also, P„ = ±{l,i}{I,X,Y,Z}® n denotes the 
Pauli group, where the number of non-trivial terms in 
the tensor product is the weight of a given E 6 P n . We 
say that a space V is stabilized by a measurement op- 
erator M with all eigenvalues A — ±1 (this includes all 
Hermitian operators in P„) if M = for any state 
\ijj) in V . We will also use the term anti- stabilized if 
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M \ip) = — \ip). A space is stabilized by a set M. of mea- 
surement operators if it is simultaneously stabilized by 
all operators in M.. A maximal space stabilized by M 
is called the stabilized space V(A4), and M. is called a 
stabilizer of V(M). The corresponding projector is de- 
noted Pm- The projector 1 — Pm corresponds to the 
orthogonal complement V {M). For a single measure- 
ment operator, M = 2Pm — 1- 

A general quantum code {{n,K,d)) is a subspace Q £ 
Hf n of dimension K, such that any detectable error ci- 
ther takes any non-zero state G Q into a state outside 
of Q, E \ip) g Q, or acts trivially on Q, E \ip) = C E |V) 
with Ce independent of \ip). A combination E\E2 of any 
two errors from a set £ of correctable errors is detectable. 
The errors are in the same degeneracy class iff e\e 2 acts 
trivially on Q. For a distance-d code, all Pauli errors of 
weight up to (d — 1) are detectable, and all Pauli errors 
of weight up to t = [(d — 1) / 2J are corr ectable§, §. 

A stabilizer code [[n, k, d]} is defined as the stabilized 
space of an Abelian group § = (G\, . . . , G n -k) of size 
2"~ fc , — 1 g" S, generated by Hermitian Pauli operators 
Gi, i = 1, . . .n — k. Explicitly, 



CWS code Q whose word operators contain a group, 



g = {|V>>:S|^> = |V>, VSeS}. 



(3) 



The logical operators Xi, Zi, i = 1, . . . , k commute with 
the code stabilizer S; they obey the usual Pauli com- 
mutation relations. These operators, along with Gi 6 § 
and the trivial il, serve as generators of the code nor- 
malizer N, a group of operators U £ P n that preserve 
the stabilizer § under conjugation, USU^ = S, U € N, 
S G §. Each correctable error E e P„ acting non-trivially 
on the code anti-commutes with at least one generator 
Gi, and correctable errors in different degeneracy classes 
anti-commute with different subsets of §. The corrupted 
code E(Q) = {E\ip) : \ip) e Q} is anti-stabilized by those 
generators Gi that anti-commute with E. Thus, a stabi- 
lizer code can be corrected by measuring the generators 
d; the corresponding set of eigenvalues A, = ±1 forms 
the syndrome of the error. 

A codeword-stabilized (CWS) code|| f?J ((n,K,d)) is 
denned in terms of a stabilizer state \s) (which is an 
[[n, 0]] stabilizer code), and a set of K mutually commut- 
ing codeword operators W = {W^l^j C P„. Explicitly, 



Q = span({K)}f =1 , \ Wi ) = W l \s)) 



(4) 



The stabilizer S = (Si, . . . , S n ) of the state \s) is the 
maximal Abelian subgroup of the Pauli group such that 
— 1 ^ §; in the context of CWS codes it is called word 
stabilizer^. 

A CWS code is a stabilizer code iff the K word oper- 
ators Wi form an (Abelian) groupQ|. Such a CWS code 
is called additive; in this case K = 2 k with integer k. 

A union stabilizer (USt) codeB can be defined as a 



a, 6 {0,1}}. (5) 



i=l 



Here gi are generators of the group W = (gi,---,gk) 
forming an additive code Qq = span ({M / |s)}vkgw) with 
dimension Kq = 2 . The operators tj form a set T of 
m translations for the code Qq. The translated spaces 
tj(Qo) are mutually orthogonal, which implies that the 
dimension of the code Q is K — m 2 k . 

The standard form of a CWS code|| Q| is denned in 
terms of a graph Q with n vertices and a classical code C 
containing K binary codewords Ci of length n. The graph 
adjacency matrix R £ {0,1}™ X ™ defines the generators 
of the stabilizer, S t = XiZ^ 1 Z^ 2 . . . , while the 
classical codewords define the codeword operators Wi = 
Z Ci = Z\ ,x . . . Z^ n . Most importantly, the graph relates 
the error-correction properties^ of the quantum CWS 
code Q = (G,C) and the classical code C. Indeed, the 
action of a singlc-qubit error Xi on the code is equivalent 
(up to an overall phase) to that of XiSi = Z^ 41 . . . Z^ in . 
Any Pauli operator E = Z V X U can thus be mapped (up 
to a phase) to the operator Z cie ^ . Here the function 



Clg(E) = v + u-lKi + u 2 R 2 



u n Hn (mod 2) (6) 



defines the graph-induced classical (binary) map of E. 
This function also defines the error degeneracy classes of 
the code Q: two correctable quantum errors E\ and E2 
are mutually degenerate iff Clg^i) = G\g(E 2 ) 0. 
For a pair of correctable errors E±, E2 from different 
degeneracy classes, C\g(Ei) ^ Clg(i?2), the corrupted 



spaces are always orthogonal, Ei(Q) L E 2 {Q) [12 



Any CWS code is locally Clifford-equivalent to a code 
in standard form[0. For CWS codes in standard form, 
we will denote the corresponding set of word operators 
and word stabilizer as Wg and Sg, respectively. 

Exhaustive screening for CWS codes. We can 
detect errors by measuring the operator Mq = 2Pq — 1, 



p Q = in^i^ 1 - 

wew 



(7) 



The corresponding ancilla measurement circuit which 
uses 2K[n 2 + 0(n)] two-qubit gates can be constructed 
as the special case of Eq. ([l3j) below. A different cir- 
cuit which requires up to n 2 + KO(n) two-qubit gates is 
constructed in Ref. |lj. 

The operators EMqE^ stabilize the spaces E(Q). For 
a CWS code Q, these spaces are orthogonal for mutually 
non-degenerate correctable errors E. This implies that 
an error can be located by measuring such operators for E 
from different degeneracy classes. For a t-error correcting 
code we can exhaustively test all correctable errors using 
up to B(n, t) measurements [Eq. (111)]. This bound is tight 
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FIG. 1: Measurement Mi A Mo. We use the projectors Pi = Pm ( , Qi s 1 - -Pi, 1 = 0,1, and assume PiPo = fofi- Circuit 
returns |/) = |1) PiPo \tp) + 10) (1 — Pi-Po) \ip) , which is equivalent to |/) = |1) Pa^aa/o + 10) (1 — Pm 1 am ). Intermediate results: 
|o) = |000>Pi |-0> + |100> Qi \b) = |000)PiP |V>) + |010) PiQo |V-> + 1 100) QiP + |H0) QiQo |V>), |c) = |001)PiP |V>) + 
|010) PiQo + 1 100) QiPq IV 1 } + 1 1 10) QiQo the last two gate groups disentangle the first two ancillas. 
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FIG. 2: Measurement for Mi © M . Notations as in Fig. [j] 
The result |/) = |1) (QiPo+PlQo) |V>) + I°) (P1P0 + Q1Q0) \ip) 
is equivalent to |1) Pm 1 ®m \ip) + |0) (1 - Pm-l®m ) \ip)- 



for non- degenerate codes where all linearly- independent 
correctable errors are mutually non-degenerate. 

Measurement algebra. To simplify error correc- 
tion, we will first decompose multi-qubit measurements 
using the algebra of projection operators tyl [ll| . A mea- 
surement M projects a state into the stabilized space 
V(M) or its orthogonal complement. In the following we 
assume that all measurement operators commute. 

In analogy with logical AND, let Mi A Mq denote the 
measurement that stabilizes V(M±)f]V(Mo). The circuit 
in Fig. [j] shows an implementation of this combination 
using logical operations on ancillas. A different circuit 
which requires only two ancillas is given in Rcf. 

An operation analogous to logical XOR is defined 
in terms of the symmetric difference of vector spaces 
AAB = span(Ar\B- L ,Bf)A ± ). We assume that there ex- 
ists an orthogonal basis common to all spaces. Then, the 
symmetric difference A A B A C A . . . is spanned by the 
basis vectors which belong to an odd number of subspaces 
A, B,C, . . .. We define the XOR of two commuting mea- 
surements, Mi © Mo, as the measurement that stabilizes 
V{Mi) A V{Mq). The corresponding circuit [Fig. || is 
based on the easy-to-check identity Mi © M = — MiM . 

Generally, the equality symbol will denote the equiva- 
lence between measurements. If Mi Mo = MoMi, then 

M = Mi A M & V{M) = V{Mi) n T(M a ), (8) 
M = Mi © M V(M) = V{Mi) A V(M ). (9) 

Decomposition of an additive CWS code. Con- 
sider an additive CWS code Qo with the set of word op- 
erators Wo = Wo = (gi, . ■ ■ , <7fc) forming a group. This 
code is a stabilizer code|?j]; it is the common stabilized 
space of the n — k generators G; of the code stabilizer So, 



Qo = n"=i fc T-'iGi)- According to Eq. (H), we also have 



Mo = M So = f\ d, 



(10) 



and can construct the corresponding measurement circuit 
by analogy with Fig. [l] using associativity. This requires 
2(n — k) controlled n-qubit Pauli operators and (n — k — 
1) thrcc-qubit Toffoli gates. Adding the corresponding 
complexities |l4| , we obtain the overall complexity of up 
to 2(n — k)(n + 3) two-qubit gates. 

This measurement can be done in the basis of the orig- 
inal CWS code. The n generators Si 6 P n of the word 
stabilizer S can be chosen Jl2[ to satisfy the orthogonality 
condition Sigj = (— l) Sij gjSi. Now, the k logical oper- 



ators of the code can be chosen as Xj = gj, 



Zj — Sji 



and the remaining generators of the orthogonalized word 
stabilizer can serve as the generators Gi = Si+k, i = 
1, . . . , n — k of the code stabilizer So. 

Decomposition of a USt code. Now consider a 
USt code Q with the set VV of word operators in the 
form Given the generators Gi of the stabilizer So 
of the additive subcode Qq, the generators of the trans- 
lated code tj(Qo) can be written as tjGiij. Then, the 
corresponding measurement operators [cf. Eq. (|l0|)1 



i—k 



Mi = i-Mot] 



(11) 



A ^4 

i=l 

The code Q is spanned by the orthogonal vector spaces 
Q = V{M Q ) = span{P(M,)}™ j , P(M Z ) _L V{M,) 



(12) 



which is equivalent to the symmetric difference Q 
P(Mi) A V(M 2 ) A ... A V{M m ). According to Eq. (|), 
this is also equivalent to the decomposition 



n—k 

.i=i 



(13) 



Since the XOR ( "ffi" ) of several measurements is imple- 
mented as concatenation [Fig. 0], it requires no overhead; 
the resulting complexity is then given by Theorem [l]. 
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Clustered measurements for CWS codes. For 
a i-error correcting CWS code Q, consider any subset 
of correctable errors, £' C £, and any correctable er- 
ror E not degenerate with those in £'. Then, the space 
£'{Q) = span Bg£ , E(Q) is orthogonal to E(Q). Fur- 
thermore, errors located on any t qubits (specified by 
the set of qubit indices A — {ii,...i t }) form a group 
of correctable errors E^ = (Xi, Z^^a- Thanks to the 
group property of E>t, for the set £' = Ea, we also 
have Q a more restrictive identity E(Qa) T Qa, where 
Qa = Ea(Q)- Thus, Qa is a quantum code which can 
detect errors E 6 £ not degenerate with those in Ea- 

Our clustered measurement technique is based on the 
observation that Qa is actually a USt code. Indeed, 
consider the original CWS code in standard form, Q — 
(G,C). The set of operators B A = {Z cx s( E ) ■ E E E A } 
forms an Abelian group of size 2 k < \Ea\ = 2 2 ' since the 
operators Z cig{ - E ^ obey the same multiplication table as 
E G E^ but are not necessarily independent. By con- 
struction, different elements of B^ are in different error 
degeneracy classes, therefore the spaces e,(Q) are mutu- 
ally orthogonal for different ei S U>a- The additional de- 
generate elements in E^ do not add to the span, therefore 
Qa = Ea(Q) = Oa(Q)- Since B A and Wg are combina- 
tions of Z-operators only, Qa is a USt code in standard 
form which uses the same stabilizer state \s) as Q, the 
Abelian group W = B>a, and the codeword operators Wg 
of the code Q as the translation set T [Eq. (||)] . 

To form the measurement Ma = Mq a that stabilizes 
the USt code Qa, we construct a set of (n — k) orthogonal 
generators Gi for the additive code Qo = B^span |s)), 
see Eq. (|To|) . The actual measurement [cf. Eq. (13)], 



M A = 



i — k 



(14) 



satisfies the complexity bound of Theorem [j]. The mea- 
surement Ma has eigenvalue 1 for all states in Qa, and 
— 1 for all states in Qa' 1 , which corresponds to all cor- 
rectable errors not degenerate with those in E^- 

To determine the error, we first perform measurements 
Mjp for all (but the last one) size-i index sets . Af- 
ter locating the covering cluster A with Abelian group 
B^ of size \TS>a\ = 2 s < 2 2t , we can find the error by 
going over all s < 2t subgroups of E^ with s — 1 gen- 
erators. Each measurement determines whether or not 
the omitted generator is a part of the error. The error 
is identified as a product of the generators present in all 
auxiliary codes that detected no errors. Overall, this re- 
quires up to N(n,t) measurements as in Eq. (^). Thus, 
for any code length n > 3, the former number of B(n,t) 
measurements [see Eq. (|IJ)] is reduced by a factor 



Some additional acceleration can be gained if the orig- 
inal CWS code is a USt code, with the set of codeword 
operators (§). In this case, for a given index set A, our 
scheme employs a bigger group W which includes the 
generators of both Ba and the original group W, and a 
smaller translation set T of size m < K. The complexity 
of a single measurement would then be reduced to 2mn 2 , 
compared to 2Kn 2 in Theorem |[ Screening of N(n, t) 
or fewer qubit clusters will locate the error. 

Note also that in the special case of stabilizer 
codes, our error-grouping technique is equivalent to the 
syndrome-based recovery jl^]. Indeed, for a stabilizer 
code Q = [[n, fc,d]], the degeneracy classes form an 
Abelian group E = (ei, . . . , e n _fe) whose 2 n ~ k elements 
are enumerated by different syndromes Jl5[ . To locate 
the error, we can go over all (n — k) USt codes E Q (Q) 
generated by the subgroups of E with one generator, e a , 
missing. Then, the code E Q (Q) is a stabilizer code that 
has to correct only one non-trivial error, £ a = (e Q ). The 
corresponding stabilizer S a has only one generator. Thus, 
error can be located by independent measurements of 
n — k Pauli operators, as we do to measure the syndrome. 

In conclusion, we constructed an accelerated clustered 
quantum error correction algorithm for a non-additive 
CWS code which uses a set of auxiliary USt codes asso- 
ciated with groups of correctable errors on size-t clusters. 
For a generic non-additive code, this reduces the number 
of error-correcting measurements approximately 3* times, 
compared to exhaustive screening of all correctable errors 
of weight t and smaller. 
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